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Dedicated  to  Professor  Dieter  Naumann  on  the  occasion  of  his  65th  birthday 


Only  three  well-characterized  stable  nitrogen  fluoride  oxide 
molecules,  FNO,  FNO2,  and  NF3O,  and  one  ion,  [NF20]+,  are 
known  in  which  the  fluorine  and  oxygen  atoms  are  directly 
bonded  to  the  same  nitrogen  atomJ^^  Herein,  we  report  the 
synthesis  and  characterization  of  [N3NFO]^,  a  novel  stable 
nitrogen  fluoride  oxide  cation  which  can  exist  in  the  form  of 
two  different  configurational  isomers. 

Our  recent  synthesis^^^  of  the  N5+  ion  [Eq.  (1)]  prompted 
us  to  explore  the  reactions  of  similar  cations  with  HN3.  A 
logical  candidate  for  this  type  of  chemistry  was  the  [NF20]^ 
ion.f^”^^  When  [NF20][SbF6]  was  combined  with  a  stoichio¬ 
metric  amount  of  HN3  in  anhydrous  HF  solution  between 
—45  °C  and  room  temperature,  one  of  the  fluorine  ligands  was 
replaced  by  an  azide  group  [Eq.  (2)] . 

[N2F]+  +  HN3  ^  N5+  +  HF  (1) 

[NF2O]  [SbF6]  +  HN3  ^  [N3NFO]  [SbF6]  +  HF  (2) 

The  use  of  excess  HN3  under  these  conditions  led  to  the 
formation  of  N5+[SbF6]“  and  N2O  [Eq.  (3)]. 

[N3NFO]  [SbFe]  +  HN3  ^  N5+  [SbFj-  +  N2O  +  HF  (3) 

The  mechanism  of  Equation  (3)  was  studied  experimen¬ 
tally  by  isotopic  substitution  and  theoretical  calculations 
and  involves  the  intermediate  formation  of  the  [NyO]^  ion.  The 
chemistry  taking  place  in  the  [NF20]^/HN3  system,  however,  is 
considerably  more  complex  than  that  given  by  Equations  (2) 
and  (3)  and  will  be  reported  in  a  separate  publication. 
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The  [N3NFO][SbF6]  salt  is  a  white  solid  which  is  stable  at 
room  temperature  and  is  highly  soluble  in  anhydrous  HF.  It 
does  not  detonate  when  either  struck  with  a  hammer  or 
scraped  with  a  metal  spatula,  and  deflagrates  when  heated  in 
an  open  flame.  It  ignites  ethanol,  reacts  violently  with  water, 
and  oxidizes  bromide  to  bromine.  Its  general  properties 
somewhat  resemble  those  of  N5^[SbF6]~.^^^  When  heated  to 
50-60  °C,  it  decomposes  mainly  to  [NO][SbF6]  with  evolution 
of  N2  and  NF3  [Eq.  (4)]. 

3  [N3NFO]  [SbFg]  ^  3  [NO]  [SbFg]  +  4  N2  +  NF3  (4) 

Reaction  (4)  was  studied  in  solution  by  nitrogen  NMR 
spectroscopy  using  [N3NFO]^  that  was  labeled  in  the  a  and 
Y  positions  with  resulting  in  the  formation  of  ^^N-labeled 
N2  and  unlabeled  [NO]^.  It  was  also  investigated  by  vacuum 
pyrolysis  of  the  solid.  The  Raman  spectrum  of  the  solid 
residue  showed  [NO][SbF6];  the  infrared  spectrum  of  the 
volatile  products,  trapped  at  —  196°C,  contained  NF3.  Reac¬ 
tion  (4)  most  likely  involves  an  intramolecular  attack  of  the 
fluorine  ligand  of  [N3NFO]^  on  the  electron-rich  N„  atom  of 
the  azido  group,  accompanied  by  extrusion  of  the  [NO]^  ion. 

The  intermediately  formed  FN3  molecule  can  decompose 
under  these  conditions  to  N2  and  NF3  [Eq.  (5)].^^^ 

3FN3^NF3+4N2  (5) 

In  the  and  ^^FNMR  and  the  Raman  spectra  of 

[N3NFO][SbF6],  two  sets  of  peaks  were  observed,  suggesting 
the  presence  of  two  stereoisomers.  This  interpretation  was 
confirmed  by  theoretical  calculations  at  the  B3LYP,  MP2,  and 
CCSD(T)  levels  of  theory  showing  two  catenated  minimum- 
energy  structures  of  very  similar  energies  (AE  <  0.6  kcal 
moE^  at  all  levels  of  theory,  with  either  one  or  the  other 
isomer  being  favored,  depending  on  the  level  of  calculation 
and  basis  set  used).  A  third  minimum-energy  structure  was 
also  found,  derived  from  a  five-membered  N4O  ring^^^  with  the 
fluorine  atom  attached  to  a  nitrogen  atom  in  the  a  position  to 
the  oxygen  atom.  However,  its  energy  was  39.4  kcal  moE^ 
higher  than  that  of  the  catenated  Z  isomer  at  the  B3LYP/6- 
311G(2df)  level  of  theory.  The  geometries,  obtained  for  the 
catenated  Z  and  E  isomers  at  the  CCSD(T)/6-311(2df)  level, 
the  highest  level  used  in  this  study,  are  shown  in  Figure  1.  At 
this  level,  the  Z  isomer  is  favored  by  0.6  kcal  moE\  Both 
structures  are  planar  (Q  symmetry)  and  differ  mainly  in  the 
relative  positions  of  the  F  and  O  atoms,  with  the  remaining 
bond  lengths  and  angles  being  very  similar.  In  the  Z  isomer, 
the  azido  and  the  fluorine  ligands  are  cis  to  each  other,  while 
in  the  E  isomer  they  are  trans  to  each  other. 
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Figure  1.  Stuctures  of  (Z)-[N3NFO]^  and  (E)-[N3NFO]+  calculated  at 
the  CCSD(T)/6-31  lG(2df)  level  of  theory  (bond  lengths  in  A,  bond 
angles  in  degrees;  N  blue,  O  red,  F  dark  blue). 


The  NMR  chemical  shifts  and  coupling  constants 
observed  for  both  isomers  in  HF  solution  at  — 30°C  are 
summarized  in  Figure  2.  As  can  be  seen,  the  isomeric  splitting 
is  very  large  for  the  fluorine  atoms  (Ad(^^F)  =  98  ppm)  and 
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Figure  2.  Multi  nuclear  NMR  parameters  of  (Z)-[N3NFO]^  and  (E)- 
[N3NF0]^  (chemical  shifts  in  ppm).  The  values  predicted  at  the 
PBEl  BPE/6-31 1  +G(2df)  level  of  theory  are  shown  in  square  brackets. 


decreases  along  the  nitrogen  chain  with  increasing  distance 
from  the  fluorine  atom  (A(3(N1)  =  22,  Ad(N„)  =  11,  Ad(Np)  = 
3,  A(3(Ny)  =  2  ppm).  The  assignments  to  the  two  isomers  were 
facilitated  by  the  fact  that  they  were  present  in  a  different 
ratio  of  1:1.45  for  the  Zand  £  isomer,  by  the  observed 
coupling  constants,  theoretical  calculations  (the  values  calcu¬ 
lated  at  the  PBElBPE/6-311  +  G(2df)  level  of  theory  are 
given  in  brackets  in  Figure  2),  and  by  analogy  to  those  known 
for  similar  covalent  azides.  For  N„  and  N^,  the  nitrogen 
chemical  shifts  result  in  smooth  trends  when  plotted  against 
the  electron-withdrawing  power  of  the  X  group  of  XN3,  while 
the  chemical  shifts  of  Np  are  randomly  scattered  around  a 
value  of  d  =  141  ±  28  ppm  (Figure  3).  The  agreement  between 
the  observed  and  calculated  trends  in  the  nitrogen  shifts  is 
good,  particularly  if  it  is  kept  in  mind  that  the  accurate 
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Figure  3.  Plot  of  the  nitrogen  chemical  shifts  of  (Z)-[N3NFO]^,  (E)- 
[N3NF0]^  and  related  XN3-type  covalent  azides  arranged  according  to 
increasing  electronegativity  of  the  X  group;  ♦  N^,  ■  Np,  A  N^,  O  N^. 


calculation  of  chemical  shifts  for  nitrogen  atoms  containing 
free  valence  electron  pairs  is  difficult,  unless  extensive 
electron-correlation  treatment,  such  as  that  provided  by  the 
CCSD(T)  method,  is  used.^^^  The  fact  that  in  HF  solution  at 
— 30  °C  the  E  isomer  is  slightly  favored  over  the  Z  isomer,  is  in 
contrast  to  our  highest-level  prediction  for  the  free  gaseous 
ion  at  0  K  and  can  be  attributed  to  either  shortcomings  of  the 
computations  and/or  solvation  or  temperature  effects. 

The  presence  of  an  NFO  group  in  each  isomer  is 
established  by:  1)  the  ^^FNMR  spectrum,  which  shows  two 
triplets  of  equal  intensity  with  V(^^F,^'^N)  constants  of  166.2 
and  168.4  Hz  and  chemical  shifts  of  d  =  274.7  and  373.0  ppm, 
respectively,  which  are  reasonable  when  compared  to  the 
[NF20]^  signal,  observed  in  this  study  under  the  same 
conditions  at  d  =  325.2  ppm  with  V(^X^^F)  =  259  Hz ;  and 
2)  the  NMR  spectrum  which  shows  two  broad  doublets  at 
d  =  -59.5  and  -81.1  ppm,  respectively,  with  V(^'^N,^^F) 
167  Hz,  resulting  from  coupling  with  one  fluorine  atom. 
These  values  are  again  in  good  agreement  with  the  nitrogen 
shift  of  d  =  -97.7  ppm  observed  for  the  [NF20]^  ion. 

Of  the  three  nitrogen  atoms  of  the  azido  group,  only  the 
Np  atoms  were  observable  in  the  NMR  spectra  at  d  = 
—167.4  and  —170.6  ppm  for  the  E  isomer  and  the  Z  isomer, 
respectively.  For  the  observation  of  N„  (d  =  —127.8  ppm  for 
the  E  isomer,  6  =  —138.9  ppm  for  the  Z  isomer)  and  N^  (d  = 
—80.1  ppm  for  the  E  isomer,  6  =  —82.6  ppm  for  the  Z  isomer), 
^^N  substitution  was  required.  The  observed  nitrogen  chem¬ 
ical  shifts  are  consistent  with  the  presence  of  a  covalent  azido 
group  attached  to  a  highly  electronegative  ligand  (see 
Figure  3)  and  a  nitrogen  atom  attached  to  a  fluorine  atom 
and  an  oxygen  atom. 

The  second  method  of  identification  used  for  [N3NFO]^ 
was  vibrational  spectroscopy.  The  observed  Raman  and 
infrared  spectra  (see  Tables  1  and  2)  are  in  good  agreement 
with  the  spectra  calculated  for  the  two  isomers  at  the  B3LYP, 
MP2,  and  CCSD(T)  levels  of  theory.  As  expected,  the 
CCSD(T)  calculations  resulted  in  the  best  fit  with  the 
observed  spectra  and,  therefore,  only  these  data  have  been 
included  in  the  tables. 

As  one  might  expect,  the  frequencies  of  the  bands 
characteristic  of  the  azido  groups  are  similar  to  those 
observed  in  other  covalent  azides, and  are  also  similar 
for  the  two  isomers.  However,  three  of  the  modes  (v2,  V5,  and 
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Table  V.  Observed  and  unsealed  calculated  vibrational  frequencies  [cm  and  intensities  for  (Z)- 
[N3NF0]+'"’ 


Obsd  frequency  [rel  int] 

Calcd  frequency  (IR)[Ra]l'’l 

Mode 

Ra,  HF  soln, 

Ra,  solid. 

IR,W  solid, 

Av 

CCSD(T)/6-31G(d) 

Av 

20  °C 

20  °C 

20  °C 

'“'nM 

14|sj/15|,j[d,e] 

a' 

2253  [10] 

2252  [10] 

32 

2246  (27)  [114] 

33.2 

V2 

1771  [1.1] 

1771  [0.8] 

1765  s 

0 

1810  (387)  [6.2] 

0 

V3 

1159  [1.7] 

1 1 62  [1 .8] 

1156  vs 

12-20'^] 

1 201  (285)  [49] 

20.5 

V4 

979  [1] 

980  sh,  vw 

994  ms 

991  (71)  [0.40] 

13.3 

V5 

733  [0.8] 

732  m 

750  (58)  [8.3] 

4.4 

Ve 

680  sh 

683  [1] 

_[f| 

669  (1.7)  [9.3] 

0.7 

V7 

546  [1.3] 

548  [0.6] 

545  w 

538  (6.7)  [5.6] 

1.3 

Vg 

443  [0.6] 

447  [0  +  ] 

436  (0.7)  [0.87] 

0.6 

V9 

204  [4.4] 

210  [1.5] 

191  (1.1)  [2.9] 

2.4 

a"  Vio 

_[k] 

683  [1] 

_[f| 

656  (9.3)  [0.77] 

3.0 

Vn 

461  [0.4] 

465  [0.3] 

451  (3.8)  [0.59] 

0 

V12 

140  sh 

140  sh 

142(0.37)  [0.24] 

0.8 

[a]  The  spectra  were  recorded  for  a  mixture  of  the  Z  and  E  isomers;  therefore,  some  of  the  observed 
relative  intensities  contain  contributions  from  the  other  isomer,  whenever  their  frequencies  coincide;  the 
Raman  spectra  were  recorded  over  a  temperature  range  of  —110  to  20 °C;  because  the  observed 
frequencies  changed  within  the  temperature  range  by  less  than  2  cm“^  for  all  bands,  only  the  data  at 
20 °C  are  listed;  in  addition  to  the  bands  from  the  cations,  the  following  bands  were  observed  for  [SbFg]”: 
Ra  (HF):  v^,  654,  vs;  V2,  571,  w;  V5,  284,  m;  Ra  (solid):  658,  567,  285  cm’^  IR  (solid):  V3,  651  cm-\ 

[b]  Calculated  IR  (in  parenthesis)  and  Raman  intensities  (in  square  brackets)  are  given  in  kmmol”^  and 
amu“\  respectively;  the  Raman  intensities  were  calculated  at  the  MP2/cc-pvtz  level,  [c]  The  IR 

spectrum  was  recorded  as  a  pressed  AgCl  disk,  [d]  Frequency  differences  between  the  isotopomer 
containing  in  the  N„  position  and  that  containing  in  the  position,  [e]  The  calculated  isotopic 
shifts  were  obtained  at  the  B3  LYP/aug-cc-pvdz  level  and  scaled  by  the  same  factor  by  which  the 
calculated  frequencies  deviated  from  the  observed  ones,  [f]  A  strong  band  was  observed  at  2241  cm“^ 
with  a  shoulder  on  its  high-frequency  side;  decomposition  studies  showed  that  this  band  did  not  decay 
together  with  the  other  [N3NFO]^  bands  and,  therefore,  is  attributed  to  an  unknown  impurity;  the  high- 
frequency  shoulder  probably  represents  of  the  cation,  [g]  Based  on  the  line  width  of  a  poorly  resolved 
band,  [h]  Obscured  by  a  strong  line  from  the  teflon-FEP  sample  tube,  [i]  Obscured  by  V3  of  [SbFg]”. 
[k]  Obscured  by  V4  of  [SbFg]”. 


Table  2:  Observed  and  unsealed  calculated  vibrational  frequencies  [cm  and  intensities  for  (E)- 
[N3NF0]+.'"’ 

Obsd  frequency  [rel  int]  Calcd  frequency  (IR)[Ra]f^' 


Mode 

Ra,  HF  soln, 
20  °C 

Ra,  solid, 

20  °C 

IR.'^'  solid, 
20  °C 

Av 

hn/isnM 

CCSD(T)/6-31G(d) 

Av 

itN/isNld.'! 

a'  Vi 

2253  [10] 

2252  [10] 

32 

2249  (16)  [126] 

33.4 

V2 

1705  [0.7] 

1  704  [0.6] 

1  704  ms 

0 

1  739  (273)  [8.9] 

0.5 

V3 

1159  [1.7] 

1162  [1.8] 

1156  vs 

1 2-20[s’ 

1210  (413)  [44] 

12.0 

V4 

991  [0.9] 

995  [1.3] 

994  ms 

10 

1009  (12)  [5.6] 

9.8 

V5 

901  [0.4] 

903  m 

18.7 

885  (24)  [7.2] 

18.5 

Ve 

668  [1] 

672  [0  +  ] 

658  (14)  [6.6] 

0.8 

V7 

516  (0.7) 

519  [0.6] 

518  w 

504  (5.3)  [3.9] 

1.8 

Vg 

465  [0  +  ] 

465  [0.4] 

463  w 

451  (1.7)  [0.95] 

0.9 

V9 

180  sh 

193  sh 

177  (1.1)  [4.2] 

2.0 

a"  Vio 

_[k] 

_[0 

675  (10)  [0.86] 

1.5 

Vn 

480  [0  +  ] 

477  vw 

473  (1.5)  [0.39] 

0.4 

V12 

140  sh 

140  sh 

159  (0.00)  [0.39] 

0.8 

[a]-[k]  See  footnotes  [a]-[k]  of  Table  1 . 

V7),  which  involve  mainly  motions  of  the  FNO  group,  differ 
significantly  and,  in  the  two  isomers,  are  separated  by  47, 160, 
and  31  cm~\  respectively,  in  HF  solution  and  by  67,  174,  and 
26  cm~\  respectively,  in  the  solid  state.  Additional  support  for 
the  predicted  structures  comes  from  the  good  agreement 
between  observed  and  calculated  isotopic  shifts  in  [^^N-N-N- 
NFO]+  and  [N-N-i^N-NFO]+  (Tables  1  and  2). 


The  assignment  of  the  vibra¬ 
tional  spectra  was  tricky  because  of 
the  presence  of  two  isomers  which 
could  not  be  separated,  as  well  as 
some  starting  materials,  such  as 
[NFjOl+P’  or  [H2N3]+,t“>  and  some 
decomposition  products,  such  as 
N5+,P>  [NO]+,t“>  and  [N02]+J“>  A 
large  number  of  spectra  were 
recorded,  and  the  assignments  to 
the  two  [N3NF0]^  isomers  were 
based  on  the  observed  growth  and 
decay  patterns.  A  typical  Raman 
spectrum,  in  which  [NF20]+  and 
[N02]^  were  present  as  impurities, 
is  shown  in  Figure  4.  Whereas  the 
1407  cm~^  band  of  [N02]^  and  the 
1869  and  1158  cm"^  bands  of 
[NF20]^  presented  no  problem, 
the  900  cm-i  band  of  [NF20]+  fell 
right  on  top  of  the  relatively  weak 
V5  band  of  the  E  isomer.  The  iden¬ 
tity  and  relative  intensity  of  the 
latter  were  confirmed  by  sub¬ 
stitution  in  either  the  or 

N„  positions  of  the  azido  group, 
resulting  in  a  50:50  mix  of  the  two 
isotopomers,  [^^N-N-N-NFO]^  and 
[N-N-^^N-NF0]+.  Both  isotopom¬ 
ers  shifted  away  from  the  unsub¬ 
stituted  [NF20]^  peak  and  resulted 
in  two  well-resolved  Raman  bands 
of  equal  intensity  with  a  frequency 
separation  of  18.7  cm~^  (calculated 
value  18.5  cm~^). 

Because  the  identity  of  the 
vibrations  in  the  500  to  1200  cm“^ 
region  was  not  obvious,  normal 
coordinate  analyses  were  carried 
out  for  both  isomers.  The  results 
are  summarized  in  Tables  SI  and 
S2  in  the  Supporting  Information, 
and  show  that  most  of  the  normal 
modes  are  complex  mixtures  of  the 
corresponding  symmetry  coordi¬ 
nates.  For  example,  the  N— F 
stretching  motion  in  the  E  isomer 
is  distributed  over  several  modes, 
with  the  highest  contribution  to  a 
single  mode  (in  this  case  V4)  being 
only  29%.  Therefore,  the  assign¬ 
ment  of  characteristic  modes,  such 
as  N— F  stretching,  to  single  frequencies  and  the  comparison 
of  such  frequencies  for  the  two  isomers  is  unwarranted. 
Owing  to  the  large  size  of  the  A'  blocks  and  the  strong 
coupling,  the  determination  of  explicit  F  matrices  and  inter¬ 
nal  force  constants  were  not  feasible. 

In  conclusion,  the  new  nitrogen  fluoride  oxide  cation 
[N3NF0]+  was  prepared  from  [NF20]+  and  HN3  and  charac- 
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Figure  4.  Typical  Raman  spectrum  of  a  mixture  of  solid  (Z)-[N3NFO]- 
[SbFg]  and  (E)-[N3NFO][SbF5]  (denoted  on  the  spectrum  as  Z  and  E, 
respectively),  containing  some  [NF20][SbF6]  (A)  and  [N02][SbF5]  (B)  as 
impurities. 


terized  by  vibrational  and  multinuclear  NMR  spectroscopy 
and  electronic-structure  calculations.  The  cation  is  planar  and 
can  exist  in  the  form  of  two  stereoisomers,  depending  on 
whether  the  azido  group  and  the  fluorine  ligand  are  cis 
(Z  isomer)  or  trans  {E  isomer)  with  respect  to  each  other. 
Although  the  azido  groups  in  the  two  isomers  are  similar,  the 
bond  lengths  and  the  NMR  chemical  shifts  of  the  N-F 
bonds  in  the  two  isomers  differ  significantly,  with  the 
Z  isomer  exhibiting  the  longer  bond  (by  0.05  A). 


Experimental  Section 

Caution  !  Neat  HN3  is  highly  explosive  and  should,  whenever  possible, 
be  handled  only  in  solution.  Anhydrous  HF  can  cause  severe  burns  and 
contact  with  the  skin  must  be  avoided. 

Materials  and  apparatus;  All  reactions  were  carried  out  in  teflon- 
FEP  or  Kel-F  ampules  that  were  closed  by  stainless-steel  valves. 
Volatile  materials  were  handled  in  a  stainless-steel/teflon-FEP 
vacuum  line.^^^^  All  reaction  vessels  were  passivated  with  CIF3  prior 
to  use.  Nonvolatile  materials  were  handled  in  the  dry  nitrogen 
atmosphere  of  a  glove  box.  Raman  spectra  were  recorded  directly  in 
the  teflon  or  Kel-F  reactors  in  the  range  4000-80  cm"^  on  a  Bruker 
Equinox  55  FT-RA  spectrophotometer  by  using  a  Nd-YAG  laser  at 
1064  nm.  The  infrared  spectra  were  recorded  on  a  Midac,  M  Series, 
FT-IR  spectrometer  using  AgCl  pellets.  The  pellets  were  prepared 
inside  the  glove  box  by  using  an  Econo  minipress  (Barnes  Engineer¬ 
ing  Co).  The  starting  materials  [NF2O]  [SbFg]  and  HN3p^  were 
prepared  by  literature  methods.  HF  was  dried  by  storage  over  BiF5.f^^] 

[N3NFO][SbF6]:  In  a  typical  experiment,  [NF20][SbF6] 
(0.3  mmol)  was  added  to  a  prepassivated,  thin-walled,  5-mm  o.d. 
Kel-F  ampule,  which  was  closed  by  a  stainless-steel  valve.  On  the 
vacuum  line,  anhydrous  HF  (270  mg)  was  condensed  in  at  -196  °C 
and  the  [NF2O]  [SbF^]  was  dissolved  in  the  HF  at  room  temperature. 
The  ampule  was  cooled  back  to  -196  °C,  and  a  mixture  of  HN3 
(0.3  mmol)  and  HF  (750  mg)  was  condensed  in.  The  resulting  mixture 
was  warmed  to  temperatures  between  -45  and  -20  °C,  producing  the 
desired  [N3NFO][SbF6]  salt  in  HF  solution.  If  the  neat  solid  was 


desired,  all  volatile  material  was  removed  under  reduced  pressure  at 
-30  °C.  The  product  was  characterized  by  multinuclear  NMR  and 
vibrational  spectroscopy. 

Theoretical  methods;  The  molecular  structures,  harmonic  vibra¬ 
tional  frequencies,  and  infrared  and  Raman  vibrational  intensities 
were  calculated  by  using  second-order  perturbation  theory  (MP2,  also 
known  as  MBPT(2)[“')  or  single-  and  double-excitation  coupled 
cluster  methods, including  a  noniterative  treatment  of  connected 
triple  excitations, and  also  at  the  DFT  level  by  using  the  B3LYP 
hybrid  functional, which  included  the  VWN5  correlation  func¬ 
tional.  Hessians  (energy  second  derivatives)  were  calculated  for  the 
final  equilibrium  structures  to  verify  them  as  local  minima,  that  is,  as 
having  a  positive  definite  Hessian.  The  calculations  were  performed 
by  using  the  electronic  structure  codes  GAMESS,^^^^  Gaussian  98,^^°^ 
and  ACES  II. The  second  derivatives  were  analyzed  by  using  the 
program  BMATRIX.^^^ 
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